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PREFACE  AND  SUMMARY 


The  advent  of  artificial  satellites  makes  it  desirable  to  be 
able  to  use  measurements  of  tbe  angular  dependence  of  diffusely 
reflected  sunlight  from  a  planetary  atmosphere  to  deduce  the 
physical  properties  of  the  atmosphere.  In  this  Memorandum  it  is 
shown  that  a  combination  of  invariant  imbedding  and  quasi lineariza 
tion  may  be  used  for  this  purpose.  Some  illustrative  numerical 
results  are  provided. 
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**  .Introduetign.  In  the  first  paper  in  this  series, ^  we  presented  a 
general  formulation  of  a  significant  class  of  inverse  problems  and  out¬ 
lined  the  use  of  quasilinearization  as  a  systematic  procedure  for  ob¬ 
taining  numerical  solutions.  In  this  paper  we  wish  to  indicate  the 
application  of  these  techniques  to  some  important  questions  in  radiative 
transfer.  Our  investigations  are  aimed  at  the  complex  problem  of  de¬ 
termining  the  nature  of  planetary  atmospheres  on  the  basis  of  various 
observations  of  the  angular  dependence  of  diffusely  reflected  light, 
and  are  pertinent  to  the  design  of  experiments  in  this  area.  The  meth¬ 
ods,  however,  are  applicable  to  a  wide  class  of  physical  processes. 


Statement  of  Problem.  Consider  an  inhomogeneous,  plane-parallel, 
non-emitting  and  isotropically  scattering  atmosphere  of  finite  optical 
thickness  whose  optical  properties  depend  only  upon  t,  the  optical 
height  above  the  bottom;  see  Fig.  1. 


Incident  Flux 
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Reflected  Flux 


T 


1 


< 


FIG.  1.  THE  PHYSICAL  SITUATION 


Let  parallel  rays  of  light  of  net  flux  n  per  unit  area  normal  to  their 
direction  of  propagation  be  incident  on  the  upper  surface  in  the  direc¬ 
tion  characterized  by  where  ^  is,  as  usual,  the  cosine  of  the  angle 
measured  from  the  inward  normal  to  the  surface.  We  suppose  here  that 
the  bottom  surface  is  a  completely  absorbing  barrier.  This  is  not  an 
essential  requirement. 
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Let  us  now  employ  invariant  imbedding*^  ^  Let  r(n,|Ao,T^)  denote 
the  intensity  of  the  diffusely  reflected  light  in  the  direction  cos 
and  set  RGi.J^jTj)  -  tyr.  Then  R  satisfies  the  integro-differential 
equation 

fr  ■  -  C’*'1  +  “o1 ) R  t l<Ti>  t1  +  ?  t  R(“’'‘',Ti)  ] x 

1  0 


(1) 


x  [i+| 


with  the  initial  condition  ROi.^,0)  -  0.  The  function  X(t)  is  the 
albedo  for  single  scattering. 

The  inverse  problem  we  wish  to  consider  ia  that  of  determining  the 
nature  of  X(t)  from  measurements  of  the  reflected  flux  at  various  angles. 


3.  Analytic  Formulation, 
analytical  structure  upon 


In  order  to  do  this,  we  must  impose  some 
\(t),  which  is  to  say,  we  must  add  some  in¬ 
formation  concerning  the  hature  of  the  physical  process.  To  illustrate 
how  this  is  done,  consider  the  case  where  the  medium  consists  essen¬ 
tially  of  two  layers,  each  with  constant  albedo,  separated  by  a  thin 
zone  of  rapid  transition  from  one  value  of  the  albedo  to  the  other. 
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Let  the  albedo  have  the  form 

\(t)  -  a  +  b  tanh  10(t-c) 

so  that  X^  —  a  "  b  in  Layer  1  and  X2  —  a  +  b  i-n  Layer  2. 

In  place  of  Eq.  (2) j  let  us  use  the  discrete  equations  obtained 


(2) 


from  Gauss  quadrature, 


<2h4) 
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dR. .(t.) 


N 


‘  ■  (“I1  +  “j1 )  Rij +  [* + 1  Rik<Ti>  ^  ] * 


(3) 


r1  +  l  Vi 

(T  )  —  1 

1  ^  J 

Using  these  equations,  we  generate  "observations"  by  choosing 
a  ■  .5,  b  »  .1,  c  *  .5,  and  integrating  to  a  thickness  of  ■  1.0. 

We  also  use  N  ■  7. 

Starting  with  the  values  b^  fs  r^j(l),  we  want  to  determine  the 
quantities  a,  b,  c,  and  t^,  the  thickness.  To  do  this,  we  ask  for  the 
values  a,  b,  c  and  which  minimize  the  expression 

Sml  {W  *  bu}2.  <*> 

ij 

where  rJLj(t1)  ■  ^  »olution  of  Eq.  (3). 

The  method  we  use  is  that  of  quasilinearization,  as  outlined  in 
Ref.  1  and  presented  in  detail  in  Refs.  5-7. 


I*.  Numerical  Results.  We  carried  out  three  types  of  numerical  ex¬ 
periments  : 

a.  Determine  c,  the  altitude  of  the  interface,  given  a,  b,  and 

V 

b.  Determine  the  overall  thickness,  given  a,  b,  and  c. 

c.  Determine  a,  b,  and  c,  the  two  albedos  and  the  altitude  of 
the  interface,  given  t^. 

Following  is  a  brief  tabular  indication  of  the  results.  As  will  be 
seen,  the  success  of  the  method  depends  crucially  upon  the  initial 
approximation.  In  other  words,  we  can  only  expect  the  method  to  be 
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successful  if  we  have  at  least  some  rough  idea  of  the  nature  of  the  actual 
physical  process. 


Table  1 


SUCCESSIVE  APPROXIMATIONS  OF  c,  THE  LEVEL  OF  THE  INTERFACE 


Approximation 

Run  1 

Run  2 

Run  3 

0 

0.2 

0.8 

0.0 

1 

0.62 

0.57 

2 

0.5187 

0.5024 

No 

3 

0.500089 

0.499970 

convergence 

4 

0.499990 

0.499991 

True  Value 

0.5 

0.5 

0.5 

Table  2 


SUCCESSIVE  APPROXIMATIONS  OF  XJt  X2,  AND  c 


Approximation 

Xj-a-b 

X2“a+b 

c 

0 

0.51 

0.69 

0.4 

1 

0.4200 

0.6052 

0.5038 

2 

0.399929 

0.599995 

0.499602 

3 

0.399938 

0.599994 

0.499878 

True  Value 

0.4 

0.6 

0.5 

The  latter  table  is  based  on  knowing  the  reflected  intensity  for 
seven  angles  of  reflection  for  each  of  seven  angles  of  incidence.  The 
quasi linearization  technique  involves  integrating  124  linear  differ¬ 
ential  equations  at  each  stage  and  solving  a  system  of  linear  algebraic 
equations  of  order  three.  The  calculations  took  about  two  minutes  on  an 
IIM  7044. 

Future  studies  involve  determining  the  accuracies  in  the  intensity 
measurements  required  to  yield  specified  accuracies  in  the  estimates  of 
the  properties  of  the  medium.  Experiments  for  a  quadratic  profile 
for  the  albedo  function  will  be  reported  upon  in  the  near  future. 


5 


! 

REFERENCES 


1.  Bellman,  R.  E. ,  H.  H.  Kaglwada,  and  R.  E.  Kalaba,  On  the  Identifi¬ 

cation  of  Systems  and  the  Unscrambling  of  Data-I:  Hidden 
Periodicities ,  The  RAND  Corporation,  KM-42S5-PR,  September  1964. 

2.  Chandrasekhar,  S.,  Radiative  Transfer,  Dover  Publications,  Inc., 

New  York,  1950. 

3.  Bellman,  R.  E. ,  R.  E.  Kalaba,  and  M.  C.  Prestrud,  Invariant  Imbed¬ 

ding  and  Radiative  Transfer  in  Slabs  of  Finite  Thickness. 
American  Elsevier  Publishing  Company,  New  York,  1963. 

4.  Bellman,  R.  E. ,  H.  H.  Kaglwada,  R.  E.  Kalaba,  and  M.  C.  Prestrud, 

Invariant  Imbedding  and  Time-Dependent  Transport  Processes, 
American  Elsevier  Publishing  Company,  New  York,  1964. 

5.  Kalaba,  R.  E. ,  "On  Nonlinear  Differential  Equations,  the  Maximum 

Operation  and  Monotone  Convergence,"  J.  Math.  Mech. .  Vol.  8 
No.  4,  July  1959,  pp.  519-574. 

6.  Bellman,  R.  E.  ,  and  R.  E.  Kalaba,  Quasilinearization  and  Boundary 

Value  Problems,  American  Elsevier  Publishing  Company,  New  York, 
to  appear. 

7.  Bellman,  R.  E.  ,  H.  H.  Kagiwada,  and  R.  E.  Kalaba,  "Orbit  Determina¬ 

tion  as  a  Multi-point  Boundary  Value  Problem  and  Quasilineariza¬ 
tion,"  Proc.  Nat,  Acad.  Sci.  USA.  Vol.  48,  No.  8,  August  1962, 
pp.  1327-1329. 


\ 


i 


UNCLASSIFIED 


UNCLASSIFIED 


